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ABSTRACT 

It has been recently shown that the observed morphological properties of the 
Bullet Cluster can be accurately reproduced in hydrodynamical simulations only 
when the infall pairwise velocity Vc of the system exceeds 3000 km/s (or at least 
possibly 2500 km/s) at the pair separation of 2R^ij., where i?vir is the virial radius 
of the main cluster, and that the probability of finding such a bullet-like system 
is extremely low in the standard ACDM cosmology. We suggest here the fifth- 
force mediated by a coupled Dark Energy (cDE) as a possible velocity-enhancing 
mechanism and investigate its effect on the infall velocities of the bullet-like sys- 
tems from the CoDECS (COupled Dark Energy Cosmological Simulations) public 
database. Five different cDE models are considered: three with constant cou- 
pling and exponential potential, one with exponential coupling and exponential 
potential, and one with constant coupling and supergravity potential. For each 
model, after identifying the bullet-like systems, we determine the probability 
density distribution of their infall velocities at the pair separations of (2 — 3)-Rvir- 
Approximating each probability density distribution as a Gaussian, we calculate 
the cumulative probability of finding a bullet-like system with Vc > 3000 km/s 
or > 2500 km/s. Our results show that in all of the five cDE models the 
cumulative probabilities increase compared to the ACDM case and that in the 
model with exponential coupling P{Vc > 2500 km/s) exceeds 10~^. The physical 
interpretations and cosmological implications of our results are provided. 

Subject headings: cosmology: theory — methods:statistical — large-scale struc- 
ture of universe 



^Astronomy Program, Department of Physics and Astronomy, FPRD, Seoul National University, Seoul 
151-747, Korea; jounghun@astro.snu.ac.kr 

^Excellence Cluster Universe, Boltzmannstr. 2, D-85748 Garching, Germany ; marco.baldi@universe- 
cluster.de 

•^University Observatory, Ludwig-Maximillians University Munich, Scheinerstr. 1, D-81679 Munich, Ger- 
many 



-2- 



INTRODUCTION 



The splendid success of the standard ACDM (A-Cold Dark Matter) cosmology that 
has been witnessed for the past two decades se ems to be overshadowe d by the recent dis- 
coveries of several possible anomalies (e.g., see iPerivolaropoulosI l2008l for a review). The 



so-called "Bullet Cluster" (1E0657-56 ) is one of those observational challenges that the stan- 
dard ACDM cosmology has t o face jMastropietro fc Burkert 2008 : Lee &: Komatsu 2010 



Thompson fc Nagamind l2011t lAkahori fc Yoshikawal l201ll ) . When iTucker et al.l ( 1l995l ) first 



observed it while looking for a failed cluster, it appeared as a large cloud of hot gas. Later, 
the Chandra observation revealed that it is in fact a very rare system composed of two head- 
on colliding massive clusters sX z = 0.296 in which a bullet-like subcluster is in the middle 



Markevitch et al. 


2002. 


2004; 


Clowe et al. 


2006; 


Markevitch 


2006) 



According to the results of iMastropietro &: BurkertI (120081 ) - based on high-resolution 
hydrodynamical simulations - the pairwise infall velocity of a "bullet" satellite onto the 
main cluster at a separation of 2i?vir (where /2vir is the virial radius of the main halo) has 
to exceed 3000 km/s, or possibly at least 2500 km/s, in order to reproduce the peculiar 
morphological properties of the Bullet Cluster, such as the mass ratio between the two 
colliding clusters, the large separation between the gas and CDM dist ributions, and the 
high shock speed inferred from X-r ay temperature rn e asure ments (see also ISpringel fc Farrar 
20071 ). In the hght of this result, iLee &: Komatsul (|2010| ) investigated how probable it is 
for a "bullet-like" system to have an infall velocity as high as 3000 km /s, using the cluster 
catalogs from the large- volume MICE simulations ( ICrocce et al.ll2010[ ). They found that 
the probability of find ing a bullet-like system with > 3000 km/s in a WMAP7 universe 
( iKomatsu et al.ll201ll ) is between ~ 10~^ and ~ 10^^^, which led them to conclude that the 
existence of the Bullet Cluster is ir icompatible with the ACDM cosmology. Very recently, 
Thompson fc Nagamind (120111 ) and lAkahori fc Yoshikawa (120111) have confirmed the results 
of Lee fc Komatsu (2010) and of Mastropietro fc Burkert ( 20081 ) using the data from different 
N-body and hydrodynamical simulations, respectively. 

Both of a conservative and a radical approach to the Bullet Cluster problem have been 
recently proposed. The conservative one sought for a lower infall velocity solution in a 
ACDM cosmology under the suspicion that the results from the hydrodynamical simulations 
or th e high shock velocity infer red from the X-ray temperatures may not be trustworthy 



(e.g. 



Forero- Romero et al.l 120101 ). In contrast, the radical approach attempted to figure out 
a mechanism capable of enhancing the infall velocities of th e bullet-like systems in non- 
standard cosmologies. For example, IWyman fc Khouryl ( l2010l ) claimed that in models with 
cascading gravity the probability of a bullet-like system with > 3000 km is four orders of 
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magnitude higher than in the ACDM model (see also lMoffat fc Tothll201Cll ). Yet, their result 
was based on a pure analytical speculation without any numerical backup. 

Taking the radical direction, we examine here the possibility that the morphological 
properties and the required high infall velocity of the Bullet Cluster can be attributed to the 
presence of coupled Dark Energy (cDE). In cDE models, the Dark Energy is a scalar field 
(j) with pot ential U((J)), which interacts with the CDM particles obeying the following two 
equations JWetterichlllQQsl : lAmendolalboOol . booj ): 



(j) + 3H(I) + 



dU 



Pc + 3if Pc 



(1) 



where an overdot represents a derivative with respect to the cosmic time t, H = a/a is the 
Hubble function, Mpi = l/y/SirG is the reduced Planck mass, and pcdm is the CDM density. 
The interaction between the DE scalar field and the CDM fluid determines an exchange 
of energy- momentum, and a consequent time variation of the CDM particle mass mcoM- In 
Equations ([T]) the strength of the coupling is fully determined by the coupling function 
defined as = —dlnmcDM/dcp. 

An enhancement of the infall velocity of the Bullet Cluster is naturally expected in 
cDE models since a long-range fifth-force generated by the coupling between and CDM 
particles has been found to play the role of accelerating structure formation processes (e.g., 



Mangano et al 



2008 



2003 



Baldi et al 



Maccio et aPbooi iMainini fc Bonomettdl2006l : iPettorino fc Baccigalupil 



2010 



Wintergerst fc Pettorino 201ol . and references therein). The level of 



the enhancement of the infall velocity of the Bullet Cluster will however significantly de- 
pend on the shape of the scalar self-interaction potential f/(0) as we ll as the coupling func- 
tion 0((t)). We refer to the w; i de literature on the subject (see e.g. , lAmendolal l2000l . 12004 



Pettorino fc Baccigalupil l2008l : iBaldi et al.l l2010l . l201ll : iBaldil l2011aL and references therein) 
for a more thorough description of the main basic features of cDE cosmologies. In particular, 
the specific m odel s discussed in the present paper have been fully defined and presented by 
BaldJ J2011ch andlBaldi J2011dh . 



The goal of this Paper is to study comprehensively the effect of cDE models on the 
infall ve locities of the bullet-like systems analyzing the public cluster catalogs of the CoDECS 
project (jBaldill2011cy ) - the largest suite of cosmological N-body simulations for cDE models 
to date - which includes various choices for the scalar self interaction potential U (0) and the 
coupling function /3(0). The Paper is organized as follows. In §2, we provide a brief overview 
of the N-body simulations for the cDE models and describe how to select the bullet-like 
systems from the halo catalogs. In §3, we determine the probability density distribution of 
the bullet-like systems and calculate the probability of finding a bullet cluster with infall 
velocity larger than 3000 km/s (and 2500 km/s as well) for each cDE model. In §4, we 
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provide the physical interpretation of our resuhs and discuss its cosmological imphcations. 



DATA AND ANALYSIS 



2.1. A Brief Overview of the CoDECS 



The CoDECS (Coupled Dark Energy Cosmological Simulations) project is the largest 
suite of N-body simulations for cDE models ever performed, and all its numerical ou tputs, 
including halo and sub-halo catalogs, have been recently mad e publicly available ( iBaldi 
2011dl ). Using the specific modifie d version bvlBaldi et al.l ( 120101 ) of the widely used parallel 
Tree-PM N-body code GADGET JSpringell 120051 ) . the different runs of the L-CoDECS suite 
follow the evolution of 1024^ CDM and 1024^ baryonic particles in a periodic cosmological 
box of linear comoving size 1 h'^Gpc from Zi = 99 to z = for different realizations of the 
cDE scenario. The individual CDM and baryonic particles have a mass of 5.84 x 10^° h~^MQ 
and 1.17 x lO^°/i~^M0, respectively, at ^ = 0, and the gravitational softening was set at a 
twenty-fifth of the mean linear interparticle spacing, eg = 20 h~^kpc. 

The CoDECS suite includes - besides the standard "fiducial" ACDM scenario - five dif- 



ferent cDE_cosnw]^gies: four models with an e xponential self-interaction potential (iLucchin fc Matarrese 
19851 : iRatra fc Pee"b3ll988l : IWetterichlEoSsh 



f/(0) = Ae""^ 



and one model with a SUGRA potential (IBrax &: Martini Il999l ) 

f/(0) = A^-^e"^'/' , 



(2) 



(3) 



where for simplicity the scalar field has been redefined in units of the reduced Planck mass 
Mpi. All the former models have the same value for the potential slope a = 0.08 but differ 
from one another in the coupling function: three models (named EXPOOl, EXP002, and 
EXP003) are characterized by a constant coupling /3(0) = const, with values 0.05 , 0.1 , and 
0.15, resp ectively, consistent with present observational b ounds on the DE-CDM interaction 



[see e.g. 



Bean et al.l l2008l : IXial l2009l : iBaldi fc Viell l2010l ) , while the remaining one (named 
EXP008e3) has an exponential coupling of the form /3(0) = 0.4 ■ e^*^. For the latter, instead, 
a negative constant coupling (3 = —0.15 is assumed, allowing for a non- s tandar d evolution of 



the resulting cDE model as proposed and explained in detail by Baldi (j2011c ). The combi- 
nation of a SUGRA potential and of a negative couphng determines a peculiar dynamics of 
the scalar field, that changes its direction of motion at some intermediate redshift between 
^CMB ~ 1100 and the present time. As a consequence of this inversion of motion, which for 
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the specific model presented here (named SUGRA003) occurs at Zi^v ~ 6.8, the DE equation 
of state parameter shows a "bounce" on the cosmological constant barrier = — 1. For 
this re ason this new class of cDE models has been called the "Bouncing cDE scenario" ( iBaldi 



2011d). 



All the models included in the CoDECS suit e have been normalized according to the 
latest results from WMAP7 (iKomatsu et al.ll201l[ l both for what concerns their background 
cosmological parameters at ^ = and the amplitude of linear density perturbations at 
the last scattering surface zqmb ~ 1100. The CoDECS runs are therefore fully consistent 
with present bounds on the perturbations amplitude at zqmb and allow to investigate the 
effects that a specific cDE scenario imprints on structure formation processes from zqmb to 
the present. Table [T] lists the six cosmological models included in the CoDECS suite, with 
the corresponding coupling function, potential type, and the expected value of erg that are 
derived for each scenario with the help of the linear perturbation theory. For a more de tailed 
introduction to the CoDECS project and its related public database, see iBaldil (l2011dl ). 

It is worth mentioning here that all the cDE models considered here, with the only 
exception of the SUGRA003 model, could be directly constrained by local measurements of 
the linear perturbations amplitude as at low redshifts. In particular, the EXP003 model 
might be already in tension with presently available constraints on the power spectrum 
amplitude, 0.78 < as < 0.86, which have been determined by several independent low- 
z observations such as galaxy-galaxy cor relation function, cosmic shear statistics, X-ray 



cluster abundanc e and Ly-a fores t (e.g., iMcDonald et al.l l2005l : iHetterscheidt et al.l 12007 



Henry et al.ll2009l : IWen et al.ll2010l . and references therein). On the other hand, the EXP002 



and EXP008e3 models appear still marginally consistent with this current limit on erg. 



2.2. Selecting the Bullet-Like Systems 



Using the CoDECS public halo catalogs at 2; = 0, we first construct a mass-limited 
(M > 10^^ Mq/Zi) sample of the cluster-sized halos identified in the simulations of each 
cosmological model. The Co DECS halo catalo gs have been produced by means of a Friends- 
of- Friends (FoF) algorithm (jPavis et al.lll985l ) with linking lengt h X = b x d, w here d is the 
mean inter-particle separation and the linking parameter b = 0.2 (IBaldil l2011dl ). Reapplying 
the FoF algorithm with a linking parameter of 6 = 0.33 to the mass-limited sample from the 
CoDECS halo catalogs, we find the clusters-of-clusters for each of which the most massive 
cluster is identified as the main cluster while the others are considered as satellites, as done 



m 



Lee fc Komatsul (120101 ) 
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The bullet-like systems are selected from the identifi e d clus ters of clusters according to 



the same four criteria that were used by lLee &: Komatsul (120101 ): i) the main cluster's mass 



M/i exceeds a certain cut-off value; ii) the main cluster and at least one of its satellites are 
on the way of head-on collision (| cosa| < 0.9, where a is the angle between the velocities 
of the main cluster and the colliding satellite); Hi) the mass ratio between the satellite and 
main cluster is smaller than one fifth {Mg/Mh < 1/5); iv) the satellite is located within 
the distances of (2 — 3)-Rvir from its main cluster where i?vir is the virial radius of its main 
cluster. The bullet-like systems must satisfy the first three criteria to match the observed 
properties of the Bullet Cluster (IE0657-57) while the last criterion is used to find the infall 
pairwise velocities of the main-satellite cluster pair s before the collisions. For the detailed 



explanation of these criteria, we refer the readers to iLee &: Komatsul (120101 ) 



Regarding the cut-off value of the main cluster's mass, it is worth mentioning here that 
the main cluster in the observed bullet cluster (IE0657-56) is a very massive one with mass 
of Mh > 10^^ h~^MQ. In the current CoDECS samples, however, there are too few (less than 
50) bullet-like systems satisfying Mh > 10^^ h~^MQ. To avoid the poor number statistics, we 
use lower cut-off values, Mh > 0.5 x 10^^ h-'^M^ and Mh > 0.7 x 10^^ R-^Mq. We also focus 
on the z = sample for the same reason (at higher redshifts there are too few bullet-like 
systems) . 

Table [2] lists the number of clusters {Nc) in the mass-limited sample, the number of the 
identified clusters of clusters [Ngc), the number and mean pairwise infall velocity of the se- 
lected bullet-like systems {Ni^uUet and Vc) for the two different cases of the main cluster's mass 
cutoff. As can be seen, the cDE models have systematically larger values of A^c, Ngc, Nbuiiet 
than the ACDM model, which is consistent with the picture that the dark sector coupling 
leads to a faster growth of the large-scale structures. Regarding the mean infall velocities, 
all of the cDE models except of the SUGRA003 model yield higher values of Vc- The larger 
the coupling is, the higher is the mean infall velocity of the bullet-like systems. In the 
SUGRA003 model, however, the value of Vc is slightly lower than the ACDM case when 
Mh > 0.5 X 10^^ h~^MQ, which should be attributed to the negative value of the coupling 
constant (see §4 for more discussions). 



3. INFALL VELOCITIES OF THE BULLET-LIKE SYSTEMS IN cDE 

MODELS 

For each selected bullet-like system, we measure the relative pairwise velocities, Vc = 
Yh — Vs, where Yh and Vc denote the velocities of the main and satellite clusters, respec- 
tively. Then, we bin the values of log Vc and count the number of those bullet-like systems 
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with log Vc belonging to a given differential bin, [log Vc, log Vc + d log Vc], to derive the prob- 
ability dens ity distribution of log Vc for each cDE model and for the ACDM model as well, 



as done in (ILee &: Komatsul 120101 ) . 



Figure [U plots p(\ogVc) at 2; = for the case of the standard ACDM cosmology (solid 
black histogram) and for the cases of the three cDE models with constant coupling and 
exponential potential (EXPOOl, EXP002 and EXP003 as dotted cyan, dashed green and 
dot-dashed blue histograms, respectively). As can be seen, the distribution p(logK) shows 
a tendency to develop a longer high- velocity tail in the cDE models with respect to ACDM. 
More specifically, the EXP002 model with /3 = 0.1 exhibits the longest high- velocity tail 
among the three models. 

To quantify the difference among the four distributions and to calculate the cumulative 
probability of Vc, we fit p(log Vc) to a Gaussian dist ribution by adjus t ing fo r each model the 



mean and standard deviation, {u, a^), as done in iLee fc Komatsul ( 120101 ). Figure |2] plots 
the fitting results. As can be seen, the Gaussian distribution gives a reasonably good fit 
to p(log Vc) for each case. Table [3] lists the best-fit values of (z/, a) for the three models, 
showing that the Gaussian fit to the infall velocity distribution p(log V"c) has a larger mean 
and a larger standard deviation in the constant coupling cDE models than in the ACDM 
cosmology. The infall velocity distribution with the highest mean is found, as expected, for 
the case of the EXP003 model, while EXP002 has the largest standard deviation. 

Integrating the best-fit Gaussian distributions, we also compute for each model the 
cumulative probabilities, PiVc > 3000 km/s) and PiVc > 2500 km/s), which are plotted in 
Figure |3] as a function of the constant coupling /3, where ACDM corresponds to /3 = 0. As 
can be seen, the highest values of the two cumulative probabilities are found for the case of 
the EXP002 model. The values of P(K > 3000 km/s) and P(K > 2500 km/s) increase by 
a factor of 10^ and 10*^ in the EXP002 model compared to ACDM case. It is interesting to 
note that although the EXP003 model has a stronger coupling, its cumulative probabilities 
are not higher than those of the EXP002 model. 

The same type of analysis has been carried out for the other two cDE models included 
in the CoDECS suite, namely EXP008e3 and SUGRA003. The results are shown in Figure 
m where the probability density distribution p(\og Vc) computed for these models is plotted 
and and compared to the ACDM case, while Figure [5] plots the Gaussian fitting functions 
whose best-fit means and best-fit standard deviations are listed in Table [31 As can be seen, 
the probability density distribution p(log Vc) for the case of the EXP008e3 model has both 
the largest mean and standard deviation, while the SUGRA003 model yields a smaller mean 
but a larger standard deviation of p(log Vc) than ACDM. 
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The cumulative probabilities, P{Vc > 3000 km/s) and P{Vc > 2500 km/s), are plotted 
in the top and bottom panels of Figure O respectively. As can be seen, both the SUGRA003 
and EXP008e3 models have larger cumulative probabilities than ACDM. For the case of 
SUGRA003, when the lower mass limit for the main cluster is set at 0.7 x 10^^ Mq//i, the 
cumulative probabilities, P{Vc > 3000km/s) and P{Vc > 2500km/s), increase by a factor of 
10 and 10^, respectively, compared to ACDM. However, when the mass cutoff is set at 0.5 x 
10^^ Mq//i, the SUGRA003 model shows almost no enhancement. In contrast, EXP008e3 
exhibits the maximal enhancement, no matter which lower mass limit of the main cluster is 
used. The cumulative probabilities, P{Vc > 3000 km/s) and P(K > 2500 km/s), increase by 
a factor of 10^ and 10^, respectively, compared to the ACDM result. Note that in this model 
P{Vc > 2500 km/s) reaches up to above 10~^, which indicates tha t if the required infall 



veloc ity of a bullet-like system can be possibly reduced to 2500 km/s (lAkahori &: Yoshikawa 



201 ll ). then the existence of the Bullet Cluster is not such an extremely rare event in a cDE 



model like the exponential coupling scenario EXP008e3. 

Furthermore, it is worth mentioning here that the above cumulative probabilities are 
lik ely to be underestim a ted d ue to the limited volume of the CoDECS runs. As shown 



m 



Forero-Romero et al.l (120101 ). the probability of finding a bullet cluster increases almost 
linearly with the simulation box volume. In other words, the larger volume a simulation 
box has, the more probable it is to find a bullet-like system with required infall velocity. In 
fact, we find for the case of A CDM model P(Vr > 3 000 km/s) ~ 10~^^ which is an order 
of magnitude lower than what iLee fc Komatsul (120101 ) found, P{Vc > 3000 km/s) ~ 10~^^, 



from the MICE simulation of volume 27/z~^Gpc^. Given this, we suspect that the real value 
of P{Vc > 3000 km/s) and P{Vc > 2500 km/s) would be higher than our estimates 10~^ and 
10~^, respectively. 



4. DISCUSSION AND CONCLUSION 

We have carried out a systematic investigation of how the probability of finding a bullet- 
like system with the required high-infall velocity (3000 km/s or possibly at least 2500 km/s) 
changes in the context of cDE models compared to the standard ACDM cosmology, by 
analyzing the public numerical data from CoDECS, the largest N-body simulations of cDE 
models to date, with a box size of 1 comoving h~~^Gpc^. Three models with constant coupling 
and exponential potential (EXPOOl, EXP002, EXP003), one bouncing cDE model with 
constant coupling and SUGRA potential (SUGRA003), and one model with exponential 
coupling and exponential potential (EXP008e3) are considered in the CoDECS suite, besides 
the fiducial ACDM cosmology. 
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In all of the five cDE models the cumulative probabilities have been found to be higher 
than in ACDM. It has been also shown that the maximal four and five orders of magnitude 
enhancements of the cumulative probabilities are yielded by the EXP008e3 model where 
PiVc > 2500 km/s) and PiVc > 3000 km/s) reach up to above 10~^ and 10"^, respectively. 
This is a remarkable result, considering that all the CoDECS simulations share the same 
random phases in the initial conditions and the same amplitude of density perturbations at 
the last scattering surface. Furthermore, given the po ssible underestimate of the cumulative 
probabilities due to the limited volume of CoDECS ( iForero- Romero et al.ll2010l ). we argue 
that our results should be considered conservative and that therefore the existence of the 
bullet cluster should not be regarded as an extremely rare event in the EXP008e3 model. 
Yet, in the other cDE models it should be still quite unlikely to find a bullet cluster with 
the observed morphology. 

The different amplitudes of the bullet velocity enhancement in the different cDE models 
can be qualitatively understood by considering the time evolution of the linear density and 
velocity perturbations of each scenario. By taking into account only the coupling strength 
and its integrated effect on the growth of density perturbations, in fact, our finding that the 
maximal velocity enhancement is obt ained f or the exponential coupling model (EXP008e3) 
might look surprising. As shown by iBaldil ( l2011dl ) (see the right panel of Figure 2 in the 
reference), the maximal effect on the density perturbation amplitude is obtained for the 
EXP003 model that reaches at z = an enhancement of 6{z) of about 20% relative to 
ACDM, while the EXP008e3 model shows a 10% enhancement (as also described by the 
different cxg values listed in Table [T] among the models). 

However, if we also consider the velocity perturbation variable 6, defined as 



fir- + fift 



?cfic + Ob^b 



(4) 



where 9b^c = V ■ Ub,c are the individual velocity perturbation variables, Ub^c are the pecuhar 
velocities, and Qb,c the relative density parameters for baryons and CDM, respectively, we 
can separately solve for 6{z) in the different cosmologies, and compare the time evolution of 
the velocity perturbations to the ACDM case. In the linear regime, the density and velocity 
perturbations 6b,c and 6b,c for baryons and CDM in a cDE cosmology evolve according to the 
following system of coupled differential equations: 



5'. 



51 



1 



{Sw^^l^ + fir - 1) + g{l3c 



9b 



e, + - [QbSb + fic^r. 



(5) 
(6) 
(7) 
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0', 



3 



where a prime denotes a derivative with respect to the e-folding time a = In a, fi^ and Vtr 
are the DE and radiation fractional densities, respectively, is the DE equation of state 
parameter, and the three functions f{f3c), g{/3c,4>') and Tc{4>) = 1 + 4/3f(0)/3 encode the 
deviation from the standard ACDM cosmol ogy due to the D E interaction, being {f,g,r) = 



(0, 0, 1) for an uncoupled DE field (see e.g. lAmendolal |2004| . for a more detailed description 



of the perturbation equations and a definition of the functions / and g). 



By solving the system ([MH]) we can separately compute the time evolution of the indi- 
vidual density and velocity perturbations for baryons and CDM and derive the total velocity 
perturbation 6 given by Equation (jl]). In Figure [7| we plot the ratio of the velocity per- 
turbation 6 to the ACDM case for all the models considered in the present work. As one 
can see from the figure, the situation is significantly different with respect to the evolution 
of the density perturbation 6. In particular, we display here for the first time the redshift 
evolution of the velocity perturbation 6' in a significant number of cDE scenarios, and our 
results show how a given cDE model can have a significantly different relative impact on the 
density and velocity perturbations. Figure [7] clearly shows that the maximal enhancement 
of velocity perturbations with respect to ACDM is realized by the EXP008e3 model (orange, 
long-dashed line) - consistently with our findings on the infall velocity in the bullet-like 
systems identified in the CoDECS simulations - even if the corresponding enhancement of 
the density perturbations amplitude (i.e. of as) is significantly smaller than in the EXP003 
model (blue, dot-dashed line). 

Although the solution for 6 displayed in Figure [7] is strictly valid only in the linear regime 
(as for the case of the density perturbation 6), the relative impact of the different cDE models 
on the velocity perturbations can be considered reliable also in the nonlinear regime for 
situations where the velocity vector is aligned with the gradient of the gravitational potential, 
since in such a case the extra-friction and the fifth-force contributions arising as a consequence 
of the DE interaction are also aligned with each other and keep behaving as for the linear 
case. This is in particular the situation of the bullet-like systems that we are considering 
in the present work, due to the head-on collision between the two clusters. Therefore we 
can consider the results obtained for the linear solution of the velocity perturbatio n 6 as 



indicative of the impact of the coupling even in these highly nonlinear systems. See iBaldi 



(j2011bl ) for a discussion on the linear and nonlinear effects of the friction term. 



Particularly interesting is also the case of the bouncing cDE model SUGRA003: for this 
scenario, in fact, the velocity perturbation 6 follows a similar path as for the EXP003 model 
at high redshifts, but is suddenly slowed down in correspondence to the bounce of the DE 
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scalar field at Zi^^ ^ 6.8, and is the only model to have a smaller velocity perturbation 6 as 
compared to ACDM at low redshifts. This is consistent with the background dynamics of 
the bouncing cDE scenario, that features an inversion of the scalar field motion at Zinv with 
the consequent change of sign of the friction term 0/3(0), that accelerates particles in their 
direction of motion before 2;inv, while it decelerates particles after zi-a^,. The suppression of the 
velocity perturbation variable 6 after z-^-^^ is therefore an expected effect due to the pecuhar 
dynamics o the bouncing cDE scenario. This is reflected in the lower value of the mean 
pairwise infall velocity detected in our sample for the Bouncing cDE model (see Table [3]), 
unique among the cDE models considered in our analysis. 

Nevertheless, as already mentioned above, the effect of the friction term become more 
complex in the nonlinear regime, due to the relative orientation between the local velocity 
field and the gravitational potential gradient. We speculate here that this additional degree 
of freedom that modulates the efficiency of the friction term in slowing down CDM particle 
is responsible for the larger variance detected in the pairwise infall velocity distribution in 
our sample (see again Table [3]) which is reflected in a slightly enhanced detection probability 
even for the SUGRA003 scenario, in spite of the reduction of the mean infall velocity. 

To see whether or not the enhanced velocities of the bullet-like systems in cDE models 
are really related to the enhanced values of the linear velocity perturbations, we compare the 
ratio of the mean infall velocity of the bullet- like systems for each cDE model to that for the 
ACDM model, VJV^, ,ACDM (from Table [2]) with the values of 6/6acdm at z = (from Figure 
[7]). The results are plotted in Figures [S] and [H] for the cases of > 0.5 x 10^^ h~^MQ and 
Mh > 0.7 X 10^^ H^^Mq, respectively. In each Figure, the dotted line corresponds to the case 
that the two ratios have the same values. As can be seen, the enhancements of Vc relative to 
the ACDM case indeed follow the pattern of the enhancements of 6 relative to ^acdm- In the 
SUGRA003 model, however, K/K;,acdm is significantly larger than 9/9acdm at 2; = 0, which 
reflects the peculiar evolution of the SUGRA003 velocity perturbations shown in Figure [71 

To conclude, we have performed a comprehensive investigation of the effect of coupled 
dark energy models on the pairwise infall velocity of the bullet-like systems identified in the 
public catalogs of the CoDECS N-body simulations. Since our results have clearly demon- 
strated that coupled dark energy has the effect of speeding up the bullet-like systems and 
that the strength of the effect significantly depends on the shape of the scalar potential as 
well as of the coupling function, we conclude that the Bullet Cluster will become a valuable 
constraint of the coupled dark energy scenario. 
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Fig. 1. — Probability density distribution of the infall velocities, log Vc, of the bullet-cluster- 
like systems measured within 2 < r/i?2oo < 3 at ^ = for the four different dark energy 
models: ACDM, EXPOOl, EXP002, and EXP003 as dashed, dotted, solid and dot-dashed 
histograms, respectively. The main cluster masses are M^^in > 0.7 x 10^^ Mq, for each 
case. 
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Fig. 2. — Gaussian- fit (thin solid line) to the probability density distribution of the infall 
velocities, log T^, of the buUet-cluster-hke systems measured within 2 < r/R2m < 3 at z = 
(thick solid histogram) for the four different dark energy models: ACDM, EXPOOl, EXP002, 
and EXP003 in the top-left, top-right, bottom-left and bottom-right panel, respectively. The 
main cluster masses are Mmain > 0.7 x 10^^ Mq, for each case. 
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Fig. 3. — Cumulative probabilities of the infall velocities of the bullet clusters versus the 
cDE couphng constant, /3 for the ACDM, EXPOOl, EXP002, and EXP003 models. 
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Fig. 4.— Same as Figure [D but for the SUGRA003 and EXP008e3 models (dotted and 
dashed histograms, respectively) in comparison with the ACDM case (solid histogram). 



- 20 - 



10.0 



1.0 



0.1 



ACDM 




simulanon 
Goussioln fit 



_i I i_ 



SUGRA003 




EXP008e3 



_l I l_ 



2.8 3.0 3.2 3.4 3.0 3.2 3.4 3.0 3.2 3.4 
log[vy(km/s)] log[vy(knn/s)] log[vy(knn/s)] 



Fig. 5. — Same as Figure [2] but for the SUGRA003 and EXP008e3 models in comparison 
with the ACDM case. 
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Fig. 6. — Cumulative probabilities of the infall velocities of the bullet clusters versus the 
cDE couphng constant, /3 for the SUGRA003 AND EXP008e3 models. 
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Fig. 7. — Velocity perturbation 6 for all the cDE models considered in the present work 
relative to the ACDM case 6'acdm- The vertical dotted line indicates the redshift at which 
the Bouncing cDE model SUGRA003 features the inversion of the scalar field motion, corre- 
sponding to 2;inv ~ 6.8. Note that when this transition occurs at 2;inv the trend of the linear 
velocity perturbation of the SUGRA003 model with respect to ACDM case is also inverted. 
On the other hand the fast growth of the coupling function at low redshifts for the EXP008e3 
model leads to the steep increase in the velocity perturbation relative to the ACDM case. In 
consequence the EXPOOSeS model has a larger value of ^/^acdm at z = than the EXP003 
model even though the overall amplitude of the linear density perturbation is higher in the 
EXP003 model. 
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Fig. 9.— Same as Figure E] but for the case of Mt > 0.7 x 10^^ H-^Mq. 
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Table 1. Defining properties of tfie six cosmological models (see lBaldill2011dl . for the 

detailed descriptions of the models). 



Model 


f/(0) 


a 


m 


0-8 


ACDM 


const. 






0. 


809 


EXPOOl 




0.08 


0.05 





825 


EXP002 




0.08 


0.1 


0. 


875 


EXP003 




0.08 


0.15 


0. 


967 


EXP008e3 




0.08 


0.4e^'^ 


0. 


895 


SUGRA003 


0— e<^'/2 


2.15 


-0.15 


0. 


806 
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Table 2. Model, number of clusters, number of clusters of clusters, number and mean 
infall velocity of bullet-like systems for the two different cases of the main cluster masses 



model 




Nsc 


-^bullet ^ 




[km/s] 


[km/s] 


ACDM 


121936 


22019 


184 


103 


999.56 


1073.33 


EXPOOl 


124467 


22448 


195 


116 


1037.38 


1118.27 


EXP002 


131880 


23648 


231 


155 


1069.03 


1247.71 


EXP003 


141915 


25223 


357 


240 


1112.55 


1197.93 


EXP008e3 


132992 


23795 


233 


159 


1146.49 


1221.66 


SUGRA003 


127651 


23216 


188 


106 


990.53 


1077.21 



^For Mh > 0.5 X 10^^ h'^M^ 
^For Mh > 0.7 X 10^^ h'^M^ 
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Table 3. Model and corresponding best-fit values of the mean and standard deviation of 

the Gaussian fitting function 



model 






(z/,( 




ACDM 


(2.99, 


0.076) 


(3.03, 


0.066) 


EXPOOl 


(3.00, 


0.081) 


(3.04, 


0.071) 


EXP002 


(3.01, 


0.091) 


(3.06, 


0.081) 


EXP003 


(3.04, 


0.086) 


(3.08, 


0.076) 


EXP008e3 


(3.05, 


0.091) 


(3.08, 


0.091) 


SUGRA003 


(2.98, 


0.081) 


(3.02, 


0.086) 



=^For Mh > 0.5 X 10^^ h'^Mc^ 
^For Mh > 0.7 X 10^^ h'^M^ 



